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FIRE  INDUCED  FLOWS  IN  CORRIDORS  --  A REVIEW 
OF  EFFORTS  TO  MODEL  KEY  FEATURES 

K.D.  Steckler 


ABSTRACT 

A literature  review  was  undertaken  to  identify  engineering  formulas  or  models 
which  can  be  used  to  predict  key  features  of  the  corridor- filling  process.  The 
results  of  that  review  are  presented  and  assessed.  The  filling  process  is 
viewed  as  a series  of  three  events:  a forward  gravity  current  moving  away  from 
the  fire  source,  a reflected  or  return  gravity  current  moving  toward  the 
source,  followed  by  uniform  filling  of  the  entire  corridor.  Recommendations 
for  estimating  the  filling  process  during  each  of  these  stages  are  presented. 


1 . INTRODUCTION 

Corridors  often  serve  as  conduits  for  transporting  life-threatening  combustion 
products  from  localized  fires  to  other  areas  of  a building.  For  example, 
combustion  products  from  a fire  developing  in  an  office  often  pass  through  a 
doorway  or  other  penetration  to  an  adjacent  corridor  where  they  establish  a 
buoyancy- driven  flow  beneath  the  ceiling.  Figure  1 illustrates  the  basic 
features  of  this  flow  as  outlined  by  Zukoski  [1].  Initially  a "forward" 
gravity  current  of  heated  gas  propagates  away  from  the  fire  source  along  the 
ceiling  of  the  corridor.  Since  buoyancy  drives  the  flow,  any  loss  of  buoyancy 
due  to  heat  transfer  from  the  hot  gas  reduces  the  velocity  of  the  front  and 
also  increases  its  depth.  Upon  striking  the  end  wall,  the  flow  reverses  and 
forms  a "return"  current  or  wave  beneath  the  existing  layer.  Eventually  the 
waves  disappear,  and  the  layer  tends  to  fill  (thicken)  uniformly. 

This  scenario  suggests  that  the  corridor  filling  process  can  be  modeled  in  two 
phases;  a relatively  short  initial  phase  in  which  horizontal  flow  phenomena 
are  significant,  followed  by  a much  longer  uniform- filling  phase  during  which 
horizontal  effects  are  small.  The  latter  is  considerably  simplier  and  easier 
to  model.  Indeed,  the  uniform- filling  process  is  precisely  the  process 
modeled  by  well-established  control -volume  or  "zone"  fire-growth  models  [2-4], 

The  initial  phase  is  less  well  understood.  For  some  applications,  such  as 
fire -resistance  assessment,  this  is  a moot  point  because  the  initial  phase  can 
be  neglected.  In  others,  it  cannot  be  neglected  because  the  time  scale  of  the 
gravity  current  phenomena  is  not  small  relative  to  the  time  scale  of  the 
process  of  interest.  Gravity-current  fronts  are  known  to  propagate  away  from 
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large  fires  at  initial  speeds  typically  in  the  range  0.2  to  2.0  m/s1 *  [1,5]. 

For  "typical"  corridors  having  moderate  to  long  lengths,  the  time  scale  for 
the  front  to  transit  the  length  of  the  corridor  could  be  on  the  order  of 
minutes.  During  this  period,  there  could  be  dramatic  variations  in  conditions 
along  the  length  of  the  corridor  which  have  a bearing  on  fire -protection 
issues  such  as  detection  and  escape  time. 

The  objective  of  the  current  work  is  to  identify  or  develop  engineering 
formulas  or  models  which  predict  the  key  features  of  the  corridor  filling 
process.  These  predictive  methods  combine  gravity- current  models  for  the 
relatively  short  initial  stage  and  zone  models  for  the  subsequent  and  longer 
uniform- filling  stage. 


2.  FORWARD  GRAVITY  CURRENT 

Benjamin's  [7]  work  on  general  gravity  currents  serves  as  a basis,  or  at  least 
a reference  point,  for  studies  of  fire -induced  forward  gravity  currents.  He 
found  that  for  an  inviscid  fluid  without  mixing  the  front  or  "nose"  of  an 
energy- conserving  gravity  current  moved  into  a stagnant  fluid  with  a velocity, 
vf , given  by  the  expression 

vf  = (g  ApS/2Pco)*  (1) 

where  g is  gravitational  acceleration,  S the  depth  of  the  current,  pm  the 
density  of  the  stagnant  fluid,  A p=pc0-p,  and  p is  the  density  of  the  gravity 
current.  Benjamin  expanded  the  theory  to  account  for  energy  dissipation  due 
to  friction  at  the  wall  or  turbulence.  Hinkley  [8]  applied  the  latter  theory 
to  fire -induced  corridor  flows.  Using  the  ideal  gas  law  to  express  the 
density  difference  Ap  in  terms  of  the  temperature  difference  AT^-T^,  he 
obtained 

vf  = CU' { (H-6) (2H-5)/[H(H+S) ] }1/3  (2) 

where  U'  = [ gq ' T/(cp p^T*  ) ] 1 ' 3 (3) 

and  C is  an  empirical  constant  to  account  for  dissipation,  H is  the  height  of 
the  corridor,  q'  is  the  heat  flow  rate  in  the  layer  at  the  front  per  unit 
corridor  width,  T is  the  absolute  temperature  of  the  layer,  Tm  the  temperature 
of  the  stagnant  layer  having  density  p^ , and  cp  is  the  specific  heat  of  the 
stagnant  layer.  Benjamin  found  C to  be  approximately  0.8.  The  velocity  of 
the  front  can  also  be  expressed  in  terms  of  the  mass  flow  rate  to  the  layer, 
m , as 


1It  is  important  to  note  that  the  present  review  addresses  buoyancy- 
driven  non- reacting  gravity  currents  during  the  early  stages  of  a fire. 
Corridor  flow  velocities  greater  than  2 m/s  can  occur  in  fire  situations  in 
which  forced  and/or  reacting  flows  are  present  in  the  corridor.  Emmons  [6] 

addresses  some  of  these  issues  including  the  possibility  of  a rapidly  moving 

smoke  pulse  due  to  roof  collapse. 
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vf  = m/(pWS)  = mT/tp^WS) 


(4) 


where  W is  the  width  of  the  corridor.  Eqs . (2-4)  allow  the  velocity  of  the 
front  and  the  depth  of  the  layer  to  be  calculated  when  the  heat  flow  to  the 
layer  per  unit  width,  the  mass  flow  to  the  layer,  and  the  temperature  of  the 
mass  flow  are  known.  It  is  important  to  note  that  Eqs.  (2-4)  assume  constant 
q'.  constant  m.  and  no  heat  loss  from  the  layer, 

Hinkley  also  examined  the  case  in  which  heat  is  transferred  from  the  hot 
gravity  current  to  the  ceiling  and  walls  of  the  corridor.  As  mentioned  above, 
heat  transfer  will  decrease  the  frontal  velocity  as  the  front  progresses  down 
the  corridor.  Assuming  no  variation  with  height  of  the  temperature  within  the 
gas  layer,  invoking  the  Boussinesq  approximation  (p-^p^  , except  in  the  buoyancy 
term),  and  using  the  notation  displayed  in  Fig. 2,  Hinkley  found 


vf  (x) 

= Vq  exp [ -K(x  - x0)] 

(5) 

Xf-x0 

- K In  [ (v0  t/K)  + 1] 

(6) 

K 

= hW/(3mcp) 

(7) 

S(t) 

= mt/lp0W(xf-x0)] 

(8) 

where  h is  the  heat  transfer  coefficient  between  the  hot  gas  and  ceiling,  t is 
time,  x is  the  horizontal  coordinate,  and  subscripts  0 and  f denote  the 
source-end  of  the  corridor  and  the  front,  respectively.  The  attendant 
temperature  variation  with  position  along  the  length  of  the  corridor  is  given 
by 


T(x)  = Tffl  + (T0-Tffl)  exp [ - 3K(x-x0 ) ] ; x < xf  (9) 

where  T is  the  temperature  of  the  gravity  current,  T0  is  its  temperature  at 
xQ  , and  T,,,  is  the  temperature  of  the  ceiling  and  remaining  air  in  the 
corridor.  Note  that  T0 , Tm , m and  q'  are  assumed  constant.  The  velocity  vQ 
appearing  in  Eq . (5)  can  be  obtained  from  Eqs.  (2-4).  Equations  (5-9)  apply 
to  "warm"  (Boussinesq)  gas  layers  which  transfer  heat  to  the  ceiling. 

For  "hot  gas",  Hinkley' s analysis  can  be  extended  to  yield 

vf(x)  = v0  { (Tra/T0 )exp[ -3K(x-x0 ) ] 

+ ((T0-Ta,)/T0)exp[-6K(x-x0)]  }1/3  (10) 

<5(t)  = mt/ [ pQ  (T0 /Tm  )W|  ] (11) 


where 


| = (xf-x0)  + (l/3K)ln{(Ta,/T0)  + [l-(Ta,/T0)]exp[-3K(xf-x0])  (11a) 

Note  that,  as  required,  Eq . (10)  reduces  to  Eq . (5)  and  Eq . (11)  reduces  to 
Eq . (8)  as  Tq-^-T^  . Unlike  Eq . (5),  however,  Eq . (10)  cannot  be  easily 
integrated  to  produce  a closed-form  solution  for  xf(t).  Thus  Hinkley' s 
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approach  can  be  applied  to  the  hot- gas  case,  but  numerical  methods  are 
required. 

Heskestad  [5]  conducted  a series  of  full-scale  corridor- flow  experiments  in 
which  he  determined  forward  and  reverse  gravity  current  propagation 
velocities.  The  experimental  arrangement  was  basically  a room/corridor 
configuration  with  the  fire  located  in  the  room.  Both  steady  and  "t- squared" 
(heating  rate  of  source  increases  with  the  square  of  time)  fires  were  used. 
These  well - instrumented  experiments  provide  a data  base  for  assessing 
corridor- flow  theories. 

Heskestad  analyzed  his  forward- current  velocity  results  in  terms  of  a "slab" 
model  which  he  developed  as  part  of  the  study.  The  model  is  expressed  as 

vf  = (g  ApS/2p)*  (12) 

which  is  similar  but  not  equal  to  Benjamin's  (Eq.  (1))  in  which  pffl  appears  in 
the  denominator.  In  contrast  to  Hinkley,  Heskestad  defined  the  convective 
heat  flow  per  unit  corridor  width  at  the  source-end  of  the  corridor  as 

q'  = cp  pvf  5 AT  = CpVfSApT*  (13) 

Solving  Eq . (13)  for  6Ap  and  substituting  into  Eq . (12)  yields 

V,  - (g  q'/^.CpT,,)1'3  (14) 

which  can  be  expressed  in  terms  of  temperature  as 

vf  = { gq ' ( 1+AT/Tm ) / ( 2pm cp Tffl ) } 1 ^ 3 (15) 

In  analyzing  his  corridor  flow  data,  Heskestad  interpreted  AT  as  the  maximum 
temperature  rise,  ATm , in  the  flow  at  the  source-end  of  the  corridor.  He 
found  that  the  experimental  forward- front  velocities  could  be  predicted  quite 
well  with  the  expression 

vf  = 0.9  (g  q‘  (l+ATm/Tt0)/(2p<0CpT(0)}1/3  (16) 

Ratios  of  his  calculated  to  experimental  front  velocites  for  three  experiments 
conducted  in  a 2.44  m x 2.44  m x 18.6  m corridor  are  plotted  as  unfilled 
symbols  in  Fig.  3.  It  is  noteworthy  that  such  good  agreement  is  achieved  even 
though  the  effects  of  heat  transfer  to  the  walls  and  ceiling  of  the  corridor 
are  neglected  in  Eq . (16). 

Although  experimental  layer  depths  (5)  were  also  reported  by  Heskestad, 
corresponding  predictions  were  not.  Nevertheless,  layer-depths  could  be 
easily  calculated  from  the  available  data  and  Eq . (13).  Ratios  of  calculated 
to  experimental  layer  depths  are  shown  in  Fig.  4 (unfilled  symbols). 

Agreement  is  poor  near  the  source -end  of  the  corridor,  but  improves  as  the 
distance  from  the  source,  xf , increases. 

Heskestad' s data  also  provide  a limited  means  for  checking  Hinkley' s analysis. 
It  is  important  to  note  that,  owing  to  heat  losses  to  the  source  - room ' s 
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surfaces . Heskestad  reports  values  of  q'  near  the  source-end  of  the  corridor 
which  increase  with  time;  even  for  experiments  using  steady  fire  sources  (e.g; 

' max  /l ' m i n ~ 3 to  5 ) . Hinkley's  analysis,  however,  assumes  constant  q'  at 
the  head- end  of  the  corridor.  Consequently,  checks  on  vf  and  S are  limited  to 
data  taken  when  the  front  was  at  the  head- end  location;  namely,  the  5.3m 
position  in  Figs.  3 and  4.  Good  agreement  between  experiment  and  theory 
(i.e.,  Eqs.  (2-4))  is  shown  for  this  location  in  both  figures.  Indeed, 
agreement  at  the  remaining  downstream  locations  was  found  to  be  equally  good 
but  results  are  not  plotted  because  the  theory  is  not  strictly  applicable  and 
agreement  may  be  fortuitious.  Nevertheless,  further  study  of  the  curious 
agreement  at  downstream  locations  is  recommended. 

A fairly  rigorous  hydraulic  model  of  transient  ceiling  flow  is  being  developed 
by  Emmons  [9].  He  is  attacking  the  problem  in  pieces  each  of  which  can  be 
validated.  His  analysis  currently  applies  to  adiabatic  flows  and  accounts  for 
hydraulic  jumps  that  can  occur  depending  on  inlet  conditions.  The  associated 
heat- transfer  problem  has  been  solved  [10]  and  the  solution  will  be  added  to 
the  model  in  the  near  future.  Emmons  expects  to  produce  a model  in  the  form 
of  a computer  subroutine  within  the  next  year. 

3.  RETURN  GRAVITY  CURRENT 

The  return  or  reflected  fire -induced  gravity  current  has  not  been  studied  to 
the  same  extent  as  the  forward  current.  Return  waves  have  been  measured,  but 
theoretical  models  have  not  been  developed. 

Heskestad  reports  that  the  return  waves  in  his  experiments  traveled  at 
approximately  constant  depth  and  at  approximately  the  same  speed  as  the 

forward  fronts.  In  the  case  of  a steady  0.5  megawatt  fire,  he  found  the  depth 
of  the  forward  current  to  be  23  percent  of  the  height  of  the  corridor.  The 
subsequent  return  wave  increased  the  overall  depth  of  the  layer  to  73  percent 
of  the  corridor  height,  thus  leaving  little  clear  space  above  the  floor.  This 
example  illustrates  that  the  height  of  the  interface  between  the  hot  and  cool 
layers  in  a corridor  can  change  quickly. 

Zukoski  [1]  has  used  the  salt-water  analog  technique  to  study  corridor  flows 
in  reduced- scale  models.  This  technique  actually  models  adiabatic  flows  and 
therefore  foregoes  study  of  heat  transfer  effects  on  the  flow.  Nevertheless, 
it  provides  an  economical  means  for  gaining  insight  into  many  flow  phenomena. 
For  example,  Zukoski  and  Kubota  [11]  have  found  that  for  a given  set  of  inlet 
conditions  1)  the  forward  current  moves  essentially  at  constant  velocity,  2) 
the  return  wave  velocity  is  nearly  constant  at  approximately  80  percent  of 
that  of  the  forward  wave,  and  3)  the  depth  of  the  return  wave  is  nearly  equal 
to  the  depth  of  the  forward  wave;  that  is,  with  the  passing  of  the  return 
wave,  the  overall  depth  of  the  layer  is  almost  doubled.  Although  this  study 
was  conducted  on  adiabatic  gravity  currents,  findings  2)  and  3)  are  consistent 
with  Heskestad' s findings  for  full-scale  gravity  currents  with  heat  transfer 
(see  previous  paragraph) . 

Salt-water  experiments  [11]  also  revealed  that  the  entire  layer  between  the 
front  of  the  return  wave  and  the  end  wall  was  essentially  stagnant  (Fig.  5) . 
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This  indicates  that  the  return  front  is  "fed"  directly  from  the  layer  above. 

It  is  tempting  to  assume  that  the  same  phenomenum  occurs  for  gravity  currents 
with  heat  transfer.  If  this  were  true,  then,  as  a first  approximation 
(ignoring  mixing  at  the  front  of  the  return  wave),  the  temperature  profile  of 
the  return  layer  could  be  equated  with  the  temperature  of  the  layer 
immediately  above  it  (Eq.  (9)).  The  temperature  of  the  hot  layer,  both  before 
and  after  the  passage  of  the  return  front,  would  vary  in  the  horizontal 
direction  but  not  in  the  vertical  direction.  Unfortunately,  the  vertical 
temperature  profile  measured  by  Heskestad  at  the  mid- length  position  of  his 
corridor  indicates  that  the  return  wave  is  markedly  cooler  and  more  stratified 
than  the  layer  above  it.  Since  Heskestad  did  not  measure  velocities  in  the 
layer  behind  the  front,  the  presence  of  the  "stagnant  layer"  cannot  be 
confirmed  or  denied  from  his  data.  More  work  is  required  to  resolve  the 
"stagnant- layer"  issue. 

Indeed,  the  stagnant -layer  behavior  is  quite  different  from  the  "overturning" 
and  "sinuous"  flows  that  Zukoski  [1]  observed  during  salt-water  simulations  of 
a fire  plume  in  a room  (Fig.  6).  In  the  latter  experiments,  momentum- 
dominated  ceiling  jets  at  the  point  of  wall  impingement  produced  overturning 
and  a well -mixed  upper  layer,  whereas  buoyancy- dominated  jets  at  the  same 
point  produced  sinuous  flow  and  a stratified  layer.  Zukoski  argues  that  the 
existence  of  either  momentum- dominated  or  buoyancy- dominated  flow  at  the 
initial  point  of  wall  impingement  is  primarily  a function  of  room  geometry.  He 
supports  this  position  with  the  salt-water  work  of  Baines  and  Turner  [12]  and 
Tangren  et  al  [13].  Baines  and  Turner  propose  that  the  transition  takes  place 
when  the  ratio  of  lateral  room  dimension,  L,  to  plume  height,  Hp , is  between 
one  and  two  (Fig.  7).  Since  this  ratio  is  greater  than  two  for  typical 
corridors,  sinuous  flow  and  stratification  would  be  expected.  The  discovery 
of  the  stagnant  layer  in  the  corridor  flow  further  complicates  the  issue. 
Perhaps  the  stagnant  layer  is  a consequence  of  the  two-dimensional  nature  of 
the  corridor  flow,  whereas  the  other  phenomena  arise  from  the  three- 
dimensional  nature  of  the  flow  produced  by  an  axis -symmetric  plume  centered  in 
a room.  In  any  event,  the  mechanisms  which  effect  the  character  of  the  return 
wave  as  well  as  the  overall  stratification  of  the  layer  are  not  well 
understood  at  this  time. 


4.  UNIFORM  FILLING 

In  their  salt-water  experiments,  Zukoski  and  Kubota  find  that  once  the 
reflected  wave  front  reaches  the  head-end  of  the  corridor,  essentially  uniform 
filling  occurs  throughout  the  corridor.  Full-scale  room/corridor  experiments 
conducted  by  Cooper  et  al  [14]  exibited  near-uniform  filling  of  the  corridor 
during  most  of  the  test  period.  Since  the  temperatures  of  the  ceiling  and 
walls  of  the  corridor  increase  with  time,  the  temperature  of  the  gas  layer 
would  tend  to  become  uniform  during  the  uniform- filling  period.  Therefore, 
simple  zone  modeling  of  this  portion  of  the  process  appears  appropriate. 
Indeed,  Cooper  et  al  showed  that  the  decent  of  the  interface  between  the  upper 
and  lower  layers  in  their  room/corridor  could  be  predicted  reasonably  well 
with  a single-room  zone  model  known  as  ASET  [2].  In  this  case,  the 
room/corridor  combination  was  treated  as  one  large  room. 


6 


ASET  models  the  uniform- filling  process  depicted  in  Fig.  8.  The  room  is 
closed  except  for  a leakage  path  near  the  floor.  The  base  of  the  fire  is 
located  a distance  A above  the  floor  and  a distance  H beneath  the  ceiling. 

The  heat  release  rate  of  the  fire  as  a function  of  time,  Q(t),  is  presumed  to 
be  known.  The  non-dimensional  height  of  the  interface  above  the  base  of  the 
fire,  f=Zi /Lc , and  the  non-dimensional  temperature  of  the  upper  gas  layer, 

4>= Th/Ta  , are  obtained  as  a function  of  non-dimensional  time,  r=t/t  , by 
simultaneously  solving  the  two  differential  equations  [2] 


dr/d  r 


-c1q-c2q1/3r5/3 , 

-^q- 

0, 


o<r^r0 

-<s<r<o 

r— « 


(17) 


d<f>/dT 


^[ciq-(^-i)c2q1/3r5/3  ]/(r0-o . o<r<r0 

^Clq/(r0+r),  -6<r<o 


(18) 


where 

Lc , tc , and  Qc  are  a characteristic  length,  time,  and 
heat  release  rate,  respectively, 

A is  the  area  of  the  floor, 

Ar  is  the  fraction  of  Q radiated  from  the  combustion  zone, 

Ac  is  the  fraction  of  Q conducted  through  the  room's 
surfaces , 

pa  is  the  density  of  gas  in  the  lower  layer, 

Ta  is  the  temperature  of  the  lower  layer, 

Cp  is  the  specific  heat  of  the  gas  in  the  lower  layer, 
r0=H/Lc;  q=Q/Q0;  S=A/Lc 

ci  = [(l-Ac)Qotc]/(PaCpTaALc)  (19) 
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c2  = (0.210  tc/A)  {[(l-Ar)Q0gLc2]/(paCpTa)}1/3 


(20) 


ASET  also  predicts  species  concentrations  via  an  equation  similar  to  Eq . (18). 
Further  details  of  the  model  and  its  applications  can  be  found  in  [2]. 

Nelson  [15]  used  ASET  to  estimate  the  filling  of  rooms  connected  by  doorways 
to  a corridor  containing  a hot  gas  layer.  Each  room  was  treated  as  a closed 
space  (with  a small  leak  at  floor  level)  containing  a virtual  fire  source  that 
supplied  energy  at  a rate  equal  to  the  enthalpy  inflow  through  the  doorway  of 
the  room.  He  estimated  the  mass  flow  through  the  vent  using  formulas  from 
FIREFORM  [16],  which  is  a collection  of  convenient  computational  schemes  for 
fire-safety  applications.  The  same  approach  could  be  used  to  obtain  a first- 
cut  estimate  of  the  uniform- filling  of  a corridor  connected  to  a room 
containing  a fire.  A more  refined  estimate  (based  on  fewer  assumptions)  could 
be  obtained  with  a multiroom  zone  model  such  as  FAST  [3]. 


5 . RECOMMENDATIONS 

Although  the  current  understanding  of  corridor  flows  is  far  from  complete, 
several  engineering  formulas  or  models  are  available  for  estimating  key 
phenomena.  Since  there  is  some  empiricism  associated  with  each  of  the 
gravity- current  formulas,  and  since  the  experimental  underpinings  were 
obtained  in  moderate - length  corridors  (on  the  order  of  20  m) , it  should  be 
understood  that  the  formulas  may  not  be  valid  for  very  long  corridors. 

Bearing  this  in  mind  and  assuming  that  q' , m,  and  T0  at  the  source -end  of  the 
corridor  are  known,  it  is  recommended  that  estimates  of  corridor  flow  be  made 
as  follows: 


Initial  Phase  Forward  Gravity  Current 

* If  the  horizontal  temperature  profile,  T(x),  is  not  needed,  use 
Heskestad's  model  (Eqs.  (13-14))  to  obtain  vf . This  model  appears  to  predict 
vf  quite  well,  even  when  q*  . m.  and  T0  vary  with  time.  Its  shortcoming  lies 
in  the  fact  that  the  attendant  8,  which  is  defined  by  Eq . (13),  may  be 
incorrect  by  as  much  as  a factor  of  three.  However,  since  m is  presumed  to  be 
known,  another  8 can  be  calculated  from  Eq . (4).  One  must  recognize,  however, 
that  an  inconsistency  exists  between  the  8 on  the  right  side  of  Eq . (4)  and 
the  8 implicit  in  vf  on  the  left  side  of  the  same  equation.  Also,  numerical 
integration  will  be  required  to  determine  xf(t)  when  q'  varies  with  time. 

* If  T(x)  is  needed  and  q' . m.  and  TQ  are  constant,  use  the  appropriate 
form  of  Hinkley's  model.  The  model  for  "warm"  gas  with  heat  transfer  yields 
algebraic  expressions  for  the  lateral  temperature  gradient  within  the  layer 
(Eq.  (9))  as  well  as  the  position  of  the  front  (Eq.  (6)).  The  model  for  "hot" 
gas  with  heat  transfer  (i.e„,  Eqs.  (9-11))  is  a more  general  extension  of  the 
"warm"  model,  but  requires  numerical  integration  to  determine  xf(t). 
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* If  T(x)  is.  needed  and  q'  . m.  and  TQ  are  not  constant,  none  of  the  models 
presented  in  this  report  apply.  In  this  case,  use  average  or  maximum  values 
for  q',  m,  and  T0  in  conjunction  with  Eqs . (9-11)  to  obtain  average  or  worst- 
case  conditions. 


Initial  Phase  Return  Gravity  Current 

Little  quantitative  information  is  available  for  the  return  current. 
Nevertheless,  based  on  Heskestad's  and  Zukoski's  work,  use  a simple  first-cut 
model  with  the  following  features:  1)  the  front  of  the  return  current  travels 
at  80  percent  of  the  velocity  of  the  forward  front  (or  perhaps  80  percent  of 
the  average  velocity  of  the  forward  front  if  it  varies),  2)  the  depth  of  the 
return  current  is  equal  to  the  depth  of  the  forward  current  (i.e.  the  depth  of 
the  hot  layer  is  doubled  upon  passage  of  the  return  front) , and  3)  the  lateral 
temperature  profile  in  the  return  current  is  identical  to  that  of  the  layer 
above  it  (i.e.,  no  gradient  in  the  hot  layer  in  the  vertical  direction). 


Uniform- Filling  Phase 

Once  the  return  front  reaches  the  head- end  of  the  corrider,  use  a zone  model 
such  as  ASET  or  FAST  to  calculate  the  uniform- fill ing  process 
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7 . NOMENCLATURE 


A area  of  floor  of  room 

C empirical  constant  (=0.8) 

C specific  heat  of  gas  in  room's  lower  layer  (ASET) 
c[  defined  by  Eq . (19) 

c2  defined  by  Eq . (20) 

cp  specific  heat  of  air 

g gravitational  acceleration 

H height  of  corridor  or  distance  between  base  of  fire  and  ceiling  (Fig.  8) 

Hp  height  of  plume  (Fig.  7) 

h heat  transfer  coefficient 

K defined  by  Eq.  (7) 

L lateral  room  dimension 

Lc  characteristic  length  (usually  set  at  0.3  m) 
m mass  injection  rate  into  ceiling  layer 
Q heat  release  rate  of  fire 

Qc  characteristic  heat  release  rate  (usually  set  to  Q at  t = 0) 
q Q/Q0 , normalized  heat  release  rate 

q'  convective  heat  injection  rate  per  unit  width  into  ceiling  layer 
t time 

t characteristic  time  (usually  set  at  1 second) 

T temperature  of  gravity  current  or  layer 

Ta  temperature  of  room's  lower  layer  (ASET) 

Th  temperature  of  room's  hot  layer  (ASET) 

T0  temperature  of  flow  injected  into  ceiling  layer 
Tffl  ambient  temperature 

AT  T-Tffl , temperature  rise 

ATm  maximum  temperature  rise  measured  in  ceiling  layer  at  source-end  of 
corridor 

U'  characteristic  velocity  defined  by  Eq . (3) 

vQ  speed  of  flow  into  ceiling  layer  at  source-end  of  corridor  (Fig.  2) 
vf  propagation  speed  of  forward  gravity  current 
vr  propagation  speed  of  return  gravity  current 

W width  of  corridor 

x lateral  coordinate 

x0  lateral  coordinate  of  point  of  injection  of  flow  into  ceiling 

layer  (Fig.  2) 

xf  lateral  coordinate  of  forward  gravity-current  front  (Fig.  2) 

Zi  height  of  gas-layer  interface  above  base  of  fire  (Fig.  8) 

A height  of  base  of  fire  above  floor  (Fig.  8) 

5 depth  of  gravity  current  or  layer;  or  A/L  , nondimens ional  height  of 

base  of  fire  above  floor 

Zi  /Lc  , nondimens  ional  height  of  gas-layer  interface  above  base  of  fire 
f0  H/Lc , nondimens ional  distance  between  base  of  fire  and  ceiling 
Ac  fraction  of  Q conducted  through  the  room's  surfaces 

Ar  fraction  of  Q radiated  from  the  combustion  zone 

£ defined  by  Eq . (11a) 

p density  of  gravity  current  or  layer 

pa  density  of  room's  lower  layer  (ASET) 

p0  density  of  gas  injected  into  ceiling  layer  at  x0 
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Pa  density  of  ambient  air 
Ap  pm-p,  density  deficiency 

T t/tc  > nondimensional  time 

/"^a  ’ nondimensional  temperature  in  upper  gas  layer 
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Figure  1.  Basic  Features  of  Corridor  Flow  [1]. 
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Figure  2.  Parameters  for  Analysis  of  Gravity  Current  with  Heat  Transfer 
(no  mixing) . 
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Figure  3.  Comparison  of  Calculated  and  Experimental  Front  Velocities. 
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Figure  4.  Comparison  of  Calculated,  and  Experimental  Thickness  of  Forward 
Gravity  Current. 
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Figure  5.  Flow  Path  Following  Impingement  of  Gravity  Current  on  End  Wall. 
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Figure  6.  Room  Filling  [1], 
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Figure  7.  Effect  of  Aspect  Ratio  on  Filling  Pattern. 
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Figure  8.  Enclosure  Flow  Dynamics  Consistent  with  ASET  Model  [2]. 
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